Compositional (version 5.4)

Contour plot of the alpha-folded model in S^2: Contour plot of the \(\alpha\)-folded model in \(S^2\)

Description

Contour plot of the \(\alpha\)-folded model in \(S^2\).

Usage

fold.contour(m, s, p, a, n = 100, x = NULL, cont.line = FALSE)

Arguments

m

The mean vector of the folded model.

s

The covariance matrix of the folded model.

p

The probability inside the simplex of the folded model.

a

The value of a for the \(\alpha\)-transformation.

n

The number of grid points to consider over which the density is calculated.

x

This is either NULL (no data) or contains a 3 column matrix with compositional data.

cont.line

Do you want the contour lines to appear? If yes, set this TRUE.

Value

The contour plot of the folded model appears.

Details

The \(\alpha\)-transformation is applied to the compositional data and then for a grid of points within the 2-dimensional simplex the folded model's density is calculated and the contours are plotted.

References

Tsagris Michail and Stewart Connie (2020). A folded model for compositional data analysis. Australian and New Zealand Journal of Statistics, 62(2): 249-277. https://arxiv.org/pdf/1802.07330.pdf

See Also

alfa.contour, norm.contour, diri.contour, mixnorm.contour, bivt.contour, skewnorm.contour

Examples

# NOT RUN {
x <- as.matrix(iris[, 1:3])
x <- x / rowSums(x)
a <- a.est(x)$best
mod <- alpha.mle(x, a)
fold.contour(mod$mu, mod$su, mod$p, a)
# }